Synchronization reveals correlation between oscillators on networks 
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The understanding of synchronization ranging from natural to social systems has driven the 
interests of scientists from different disciplines. Here, we have investigated the synchronization 
dynamics of the Kuramoto dynamics departing from the fully synchronized regime. We have got the 
analytic expression of the dynamical correlation between pairs of oscillators that reveals the relation 
between the network dynamics and the underlying topology. Moreover, it also reveals the internal 
structure of networks that can be used as a new algorithm to detect community structures. Further, 
we have proposed a new measure about the synchronization in complex networks and scrutinize 
it in small-world and scale-free networks. Our results indicate that the more heterogeneous and 
"smaller" the network is, the more closely it would be synchronized by the collective dynamics. 
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The last decade has witnessed the rapid explosion of 
interest in complex networks, which are found in many 
fields as diverse as biology, technology and social organi- 
zations P, 0, S, 13] • The relation between structure and 
function becomes a key area in the study of complex net- 
works. In the study of networked-dynamics, the emer- 
gent synchronization of interacting oscillators generally 
has occupied a privileged position because of its rich ap- 
plications in variety of areas ran ging from Neuroscience 

to sociology H @, 0, I, 1 M, MM M M, EE 0, E3 . 

In neural systems, the brain dynamics is characterized 
by synchronization phenomena for a given topology of 
synapses. And in metabolic networks, hundreds of inter- 
connected biochemical reactions are responsible for the 
biomass and energy fluxes which are adjusted to optimize 
the robustness of synchronized behavior. Other example 
is the synchronized or coordinated behavior of commu- 
nication patterns in social organizations. The study of 
synchronization provides us with insights into the key 
issue: how the collective dynamics couple the relation- 
ships between the systematic function and the underlying 
topology. 

One of the most successful attempts to understand syn- 
chronization phenomena is due to the Kuramoto model 
(KM) [HQ, 

which is rich enough for many different 
contexts, including superconducting currents in Joseph- 
son junction arrays, emerging coherence in populations 
of chemical oscillators, and the accuracy of central cir- 
cadian pacemakers in insects and vertebrates, etc.. This 
model describes a population of N coupled phase oscil- 
lators which evolves in time according to the following 
dynamics 
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Here, 9i represents the phase of the ith oscillator, u)i the 
intrinsic frequency, a the coupling constant, Aij the ef- 



fective coupling between the oscillators, and £i(t) is the 
white noise due to the complicated environment, with 
expectation and variance 

<&(*)> = o, 

There are many attractors that each oscillator rotates at 
its natural frequency when incoherent. As the coupling 
strength increases, some units become resonant. At suffi- 
ciently strong coupling, when coherent, there is only one 
attractor of the dynamics, and all oscillators rotate at 
the same frequency, uJj, the average of u>i. We call it 
fully synchronized. 

The synchronization problem is solved mainly in the 
mean-field approach which, unfortunately, is not usually 
fulfilled in real systems. The nontrivial pattern of con- 
nectivity in complex networks brings the non mean-field 
properties and incorporates many new questions to the 
research of synchronization. Because of the elegant work 
of Pecora and Carroll [13] , the Master Stability Function 
(MSF) formalism allows us to analysis the dynamical syn- 
chronizability in terms of purely topological porperties, 
independent of details of unit dynamics. However, the 
MSF requires several constraints and has limited appli- 
cations that it only deals with the stability of the exact 
synchronized states, i.e., it is used to study infinitesimal 
deviations from the synchronization manifold of chaotic 
oscillators. However, many realistic dynamics are only 
close to or even far from the full synchronization. 

In this work we study the KM dynamics departing 
from the fully synchronized state. We get the analytic 
expression of the dynamical correlation between any two 
oscillators which reveals the relation between the net- 
work dynamics and the underlying topology. It could 
be used as a new algorithm to detect community struc- 
tures. Moreover, we propose a new measure about the 
synchronization in complex networks and scrutinize it in 
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FIG. 1: Top row: Illustration of the Santa Fe Institute col- 
laboration network. Bottom left: The analytic correlation 
matrix resulting from Eq. [Jj Bottom right: The simulated 
results of the correlation between pairs of oscillators from the 
Kuramoto model. The colors are a gradation between white 
(strong correlation) and black (weak correlation). 



Small- World and Scale- Free networks. 

When the KM dynamics is close to the attractor of 
synchronization, the phase differences are small and then 
the sine coupling function can be approximated linearly, 
so that we can analyze the linearized dynamics of the 
Kuramoto model in terms of the Laplacian matrix, 

^ = ay (ft 

j 

= -a^LijOj+bit), (2) 

3 

where Ly = J2 m a lm ~a i3 = <), , k, a,,, ay stands for 
the adjacency matrix (1 if nodes i and j are connected 
and otherwise), Sij is the Kronecker delta and the de- 
gree hi is defined as the number of nodes to which the 
node i is connected. Thus, we can rewrite the equation 



in terms of the normal modes, 
d-d 

^ = -a\ a d a + ( a (t), (3) 
at 

where d a = J2j tpajOj, (a = J2j ^aj^j, and tpaj denotes 
the ath normalized eigenvector of the Laplacian, A Q is 
the corresponding eigenvalue for a — 0, N — 1. Con- 
sidering is delta correlated, we can find easily that 
( a (t) is also delta correlated: (Ca(t)Cp(t')) = 25 a pS(t— £'). 

Without loss of generalization, we set er = 1, £ Q (0) = 0. 
Using stochastic calculus methods [2l|, we get 

MQ= / e-^^C^df, (4) 
Jo 

so that 

Jo Jo 

= f 2e- 2X ^ t - t ">dt' 
Jo 

= Ul-e~^) (5) 

Thus, the time scales of the dynamical relaxation of the 
eigenmodes are inversely proportional to the correspond- 
ing eigenvalues. And for large t 1/2 A Q , we have 

(V a (t) 2 ) = 1/Aa. (6) 

Transforming it back to the basis of 0i by substituting 
9i = ^2 a ipou'da, we find the steady-state correlation func- 
tion: 

N-l 

a=l 
N-l x 

= ^-^ailpaj, (7) 

where 9 = Yli 9i and (...) stands for the average over 
initial random phases. This analytic expression gives the 
dynamical correlation between any pair of oscillators just 
in terms of the properties of the underlying topology. Ac- 
tually, the right term of Eq. [7J gives the pseudo-inverse of 
the Laplacian matrix, which acts like a Green Function. 
It could be used to detect the community structure in 
various complex networks. 

One real world example with known community struc- 
ture is the collaboration network of the Santa Fe Insti- 
tute from ref. [22]. The top row of Fig. 1 illustrates 
the largest component of the collaboration graph, which 
consists of 118 scientists denoted by nodes. The edge 
is drawn between a pair of scientists if they coauthored 
one or more papers. Obviously, the scientists of different 
interests group in different communities. 
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FIG. 2: (color online). Illustration of the same hierarchi- 
cal structure in the eigenvalues spectra and degree series. It 
implies a progressive cascade of the synchronization pattern 
from the central core to the nodes with smaller degrees. The 
network data is used as the same as in Fig. Q] 



In the bottom row of Fig. 1, we illustrate the dynam- 
ical correlation between any pairs of nodes on which the 
KM oscillators is placed. The left one is the analytic re- 
sults from the relation Eq. [7] and the right is the results 
of numerical simulations from Eq. [TJ They accord with 
each other very well. In the two figures we can identify 
the distinct dynamics-based communities wherein nodes 
become strong correlated in groups, coherently with their 
topological structure. Moreover, smaller subcommunitics 
can be found clearly in these dynamics-based communi- 
ties which corresponds with the sketch map of the col- 
laboration network. Although some nodes in the same 
dynamics-based community are not directly connected, 
they are topologically equivalent and belong to the same 
structural level. Therefore, they have strong correlation 
on dynamics. The dynamics-based groups indicates the 
different functional modules. When applied to the spe- 
cific internal structure, which might be the fingerprint 
of different functional groups for social or biological net- 
works, the dynamical correlation may help us understand 
better the interplay between systemic structure and func- 
tion. 

In another way, we can use the mean-field approxima- 
tion to rewrite the Eq. [2] as: 



dt 



3 

= -^-6) + ^). 



Using the same stochastic calculus method, we get 



(8) 



(9) 



which indicates that the time scales of the dynamical re- 
laxation of the observed phases are inversely proportional 
to the degrees of the corresponding oscillators. The dif- 
ferent connectivities in the topology give rise to the cor- 
responding order of the eigenvalues, which are illustrated 
in Fig. 2. The hubs with the large eigenvalues induce the 
fast relaxations, i.e., the evolutional pattern of synchro- 
nization starts first at the hubs and then pervades almost 
the whole network in a hierarchical structure across the 
nodes with smaller degrees. And for large t ^> l/2fcj, we 
have got — 9i) 2 ) — 1/ki, which indicates that hubs 
are synchronized more closely by the collective dynamics. 

For a given realization of networks, we use R to denote 
the average variance of phase variables which character- 
izes the collective synchronized ability of the underlying 
network: 



2 >, 



(10) 



The smaller the value of R, the smaller the fluctuation 
of network is, i.e., the more synchronized . Substituting 
the Eq. [7] and considering the orthonormality of ip a i, we 
have 



R = 



1 N ~ 1 1 

-Y — 

Ct—1 



(11) 



Therefore, independent of details of unit dynamics, we 
can characterize the dynamical synchronizability of net- 
worked systems just in terms of the spectra of the Lapla- 
cian matrix, i.e., the purely topological porperties. It is 
the same as the spirit of the MSF. However, in contrast 
to the MSF, the new measure Eq. QT] considers the spec- 
trum in its entirety, not only the extremal eigenvalues. 
What this new measure emphasizes is on the fluctuation 
of the process of synchronization, rather than on rigorous 
bounds for the threshold of desynchronization. In the fol- 
lowing, we use this new quantity about synchronization, 
R to measure the synchronizable ability in small-world 
(SW) and scale-free (SF) networks, respectively. 

We first consider the Watts-Strogatz model of the SW 
networks 23J , which is constructed by a rewiring process 
on a regular ring graph. The probability of rewiring con- 
nection is controlled by a parameter p, by tuning which 
the obtained network possesses both short average dis- 
tance and high heterogeneity that nodes no longer have 
the same degrees, instead they follow a Poisson distribu- 
tion. We define S as the standard deviation of the degree 
distribution and D as the average network distance, aver- 
aged over all pairs of nodes. The upper layer of Fig. [3Ja) 
shows the dependence of S and D on the rewiring proba- 
bility p in the SW network. As p increases, the variance 
S increases, which implies that the degree distribution 
becomes more broad and heterogeneous. And as expect, 
the network distance D decreases as the rewiring proba- 
bility p (or the heterogeneity of the degree distribution) 
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the underlying topology. Moreover, it also reveal the in- 
ternal structure of networks that can be used as a new 
algorithm to detect community structures. Further, we 
have proposed a new measure about the synchronization 
in complex networks and scrutinize it in SW and SF net- 
works. Our results indicate that the more heterogeneous 
and "smaller" the network is, the more closely it would 
be synchronized by the collective dynamics. We hope our 
study can provide new insights into the understanding of 
the role of synchronization between the network struc- 
ture and function. We also expect it can provide new 
tools to detect community structure and to analyze the 
ubiquitous synchronization phenomenon. 
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FIG. 3: (a) The synchronizability of SW networks of size 
N = 1000, (k) — 8. The average network distance D and 
the new measure R decreases with the rewiring probability 
p in the same trend, while the standard deviation S of the 
degree distribution increases, (b) The synchronizability of SF 
networks of size N = 1000, fco = 4. The values of D and R 
increase with the scaling exponent 7, while the values of S 
are decreased. Since less R means more synchronizable, all 
the results in SW and SF networks indicate that the network 
would be more synchronizable as it becomes "smaller" and 
more heterogeneous. All plots are averaged over 100 realiza- 
tions. 



increases. In the lower layer of Fig. EJa), R is observed 
to decrease, which implies enhancement of the synchro- 
nizable ability. 

Secondly, we implement the same measures on the SF 
network. The model 24J is generated by randomly con- 
necting nodes, forcing each nodes i to have connectiv- 
ity ki > ko which follows the probability distribution 
P(k) ~ k~~< ' . Note that decreasing 7 increases S by in- 
ducing the longer tail in the connectivity distribution and 
smaller 7 gives rise to shorter D, as exhibited in the up- 
per layer of Fig. \5j[b). And in its lower layer, R is ob- 
served to increase as 7 is increased, which indicates less 
synchronized when less heterogeneous. 

The results shown in Fig. [3] imply that the synchroniz- 
able ability on the SW and SF networks is improved as 
the heterogeneity of the degree distribution is increased 
or as the average network distance is decreased. In other 
words, as the network becomes more heterogeneous and 
"smaller" , with S gets larger and D gets shorter, it would 
becomes synchronized more closely by the collective dy- 
namics. 

In summary, we have investigated the synchronization 
dynamics of the KM departing from the fully synchro- 
nized state. We have got the analytic expression of the 
dynamical correlation between pairs of oscillators that 
reveals the relation between the network dynamics and 
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